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Abstract Understanding galaxy formation is one of the most pressing issues in cosmol- 
ogy. We review the current status of galaxy formation from both an observational and a 
theoretical perspective, and summarize the prospects for future advances. 
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1 INTRODUCTION 

Numerical simulations of large-scale structure have met with great success. However these same 
simulations fail to account for several of the observed properties of galaxies. On large scales, ~ 
0.01 — 100 Mpc, the ansatz of cold, weakly interacting dark matter has led to realistic maps of the 
galaxy distribution, under the assumptions that light traces mass and that the initial conditions are pro- 
vided by the observed temperature fluctuations in the cosmic microwave background. On smaller scales, 
light no longer traces mass because of the complexity of galaxy and star formation. Baryon physics must 
be added to the simulations in order to produce realistic galaxies. It is here that the modelling is still 
inadequate. 

In this review, we will begin with the standard phenomenology of galaxy formation, then discuss 
methods and present the recent observational and modeling advances, finishing with a summary of the 
numerous outstanding issues in galaxy formation theory. 



2 PHENOMENOLOGY 
2.1 The luminosity function 

Theory provides the mass function of dark halos. Observation yields the luminosity function of galaxies, 
usually fit by a |Schechter| ( |1976| l function. Comparison of the two is at first sight disconcerting. One can 
calculate the M/L ratio for the two functions to overlap at one point, for a mass M* corresponding to 
L*. Define t coo \ = 3 / '2nkT '/ '[n 2 A(T)] and td yn — 3/ (\ / 32ttGp). For star formation to occur, cooling 
is essential, and the condition t coo i < t<x yn guarantees cooling in an inhomogeneous galactic halo where 
gas clouds collide at the virial velocity. One finds that 

Jtf* _ ^_ 7 J}P_ tcOO\ rpl + 20 
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Fig. 1 Role of feedback in modifying the galaxy luminosity function 



where a = e 2 /(he) and a g = Gm^/e 2 are the electromagnetic and gravitational fine structure con- 
stants. For a cooling function A(T) cx T@ , over the relevant temperature range (10 5 — 10 7 K), one can 
take (3 s» — 1/2 for a low metallicity plasma (Gnat & Sternberg 2007 \ . The result is that one finds a 
characteristic galactic halo mass, in terms of fundamental constants, to be of or der 1O 12 M |siikl|l977| l. 
The inferred value of the mass-to-light ratio M/L is similar to that observed for galaxies. This is a 
success for theory: dissipation provides a key ingredient in understanding the stellar masses of galaxies, 
at least for the "typical" galaxy. The characteristic galactic mass is understood by the requirement that 
cooling within a dynamical time is a necessary condition for efficient star formation (Fig. [TJ. 

However, the naive assumption that stellar mass follows halo mass, leads to too many small galax- 
ies, too many big galaxies in the nearby universe, too few massive galaxies at high redshift, and too 
many baryons within the galaxy halos. In addition there are structural problems: for example, massive 
galaxies with thin disks and/or without bulges are missing, and the concentration and cuspiness of cold 
dark matter is found to be excessive in barred galaxies and in dwarfs. The resolution to all of these 
difficulties must lie in feedback. There are various flavors of feedback that span the range of processes 
including reionization at very high redshift, supernova (SN) explosions, tidal stripping and input from 
active galactic nuclei (AGN). All of these effects no doubt have a role, but we shall see that what is 
missing is a robust theory of star formation as well as adequate numerical resolution to properly model 
the interactions between baryons, dynamics and dark matter. 



2.2 Star formation rate and efficiency 

In addressing star-forming galaxies, the problem reduces to our fundamental ignorance of star formation. 
Phenomenology is used to address this gap in our knowledge. Massive star feedback in giant molecular 
clouds, the seat of most galactic star formation, implies a star formation efficiency (SFE), defined as star 
formation rate (SFR) divided by the ratio of gas mass to dynamical or disk rotation time, of around 2%. 
This is also found to be true globally in the Milky Way (MW) disk. 

Remarkably, a similar SFE is found in nearby star-forming disk galaxies. Indeed, SFRs per unit area 
in disk galaxies, both near and far, can be described by a simple law, with SFE being the controlling 
parameter ( |Slik"l [T9971 |Elmegreen|[T997] l: 

SFR x DYNAMICAL TIME _ no 
SFE = GAS MASS °-° 2 - (1) 
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The motivation comes from the gravitational instability of cold gas-rich disks, which provides the 
scaling, although the normalization depends on feedback physics. For the global law, in terms of 
SFR and gas mass per unit area, SN regulation provides the observed efficiency of about 2% which 
fits essentially all local star-forming galaxies. One finds from simple momentum conservation that 
SFE = CT gas v CO oiTO*SN/£'sN tlal ~ 0.02. Here, w coo i is the SN-driven swept-up shell velocity at which 
approximate momentum conservation sets in and m*sN ~ 15OM is the mass formed in stars per SNII, 
in this case for a |Cha brier (2003 ) initial mass function (IMF). This is a crude estimator of the efficiency 
of SN momentum input into the interstellar medium, but it reproduces the observed global normalization 
of the star formation law. 

The fit applies not only globally but to star formation complexes in individual galaxies such as M51 



and also to starburst galaxies. The star formation law is known as the Schmidt-Kennicutt law (Kennicutt 



et al. |2007| l, and its application reveals that molecular gas is the controlling gas ingredient. In the outer 



parts of galaxies, where the molecular fraction is reduced due to the ambient UV radiation field and 
lower surface density, the SFR per unit gas mass also declines (Bigiel et al. 201 1 1. 




Fig. 2 Schmidt-Kennicutt laws on nearby (including Local Group galaxies as shaded regions) 
and distant galaxies, as well as Milky Way Giant Molecular Clouds (Krumh olz et al.||2012| l. 
The solid line is similar to equation ([TJ. 



For disk instabilities to result in cloud formation, followed by cloud agglomeration and consequent 
star formation, one also needs to maintain a cold disk by accretion of cold gas. There is ample evidence 
of a supply of cold gas, for example in the M33 group. Other spiral galaxies show extensive reservoirs of 
HI in their outer regions, for example NGC 6946 ( |Boomsma etal) |2008| l and UGC 2082 ( |Heald et al 



4 



Joseph Silk 1 ' 2 ' 3 , Gary A. Mamon 1 



|201 1| ). Recent data extends the Schmidt-Kennicutt law to z ~ 2, with a tendency for ultraluminous 
starbursts at z ~ 2 to have somewhat higher SFE ( |Genzel et al.| |2010| see Fig. [2j. 

A more refined theoretical model needs to take account of star formation in a multi-phase interstel- 
lar medium. One expects self-regulation to play a role. If the porosity in the form of SN remnant-driven 
bubbles is low, there is no venting and the pressure is enhanced, clouds are squeezed, and SN explosions 
are triggered by massive star formation. This is followed by high porosity and blow-out, and the turbu- 
lent pressure drops. Eventually halo infall replenishes the cold gas, the porosity is lowered and the cycle 
recommences. Some of this complexity can be seen in numerical simulations ( Agertz et al. 201 1\ . SNe 



provide recirculation and venting of gas into fountains, thereby reducing the SFE and prolonging the 
duration of star formation in normal disk galaxies. 



3.5 ™r 




M r absolute magnitude 



Fig. 3 Illustration of galaxy bimodality. The contours are the density of SDSS galaxies in 
color-luminosity space, after correction for selection effects (Baldr y~et al.||2004) . 



In fact, galaxy colors illustrate the bimodality of SFRs. Elliptical and lenticular galaxies are red, 
spirals are blue. This lyric does not hide a continuity in galaxy properties: most galaxies lie in either 
the Red Sequence or the Blue Cloud (see Fig. |3j. This suggests that star formation in galaxies is either 
ongoing or was quenched several Gyr ago. The small fraction of intermediate population, Green Valley 
galaxies suggests that some galaxies have experienced a recent quenching of their star formation. Seyfert 
galaxies have intermediate age stellar populations (Schawinski et al. 2007 1 see fig.Q and mostly lie in 



the Green Valley ( |Schawi nski 2012). This suggests that star formation is quenched by nuclear activity. 
2.3 Scaling relations 

The global properties of early-type galaxies are known to correlate: early work focussed on L ~ <r^ 



(Faber & Jackson 1976). The early work found a slope of 4 because of the inclusion of bright and 
faint galaxies. The modern work finds a slope of 5 for luminous galaxies (Mb <, —20.5, core-Seisic 
galaxies) and a slope of 2 for the less luminous spheroids, and has been distilled into the Fundamental 
Plane linking mass, mass-to-light ratio, and mean surface brightness at the effective radius ( |Bender et al. 
[19921 . ( |Benderetal.|[T992T ). 

Figure [5] shows a more modern version of the properties of early-type galaxies, to which are added 
globular clusters and clusters of galaxies. It is not yet understood what makes the continuity of the 
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Fig. 4 Ages of galaxies of different activity ( Schawinski et al. 2007 1 



global properties of massive systems fragment into two branches with ultra-compact dwarfs and globular 
clusters on one side and dwarf spheroidals on the other. 

2.4 Evolution of low mass galaxies 

The accepted solution for gas disruption and dispersal in intermediate mass and massive dwarfs (halo 
mass ~ 10 8 — 10 10 M©) is by SN feedback. SNe expel the remaini ng baryons in system s of halo mass 



up to ~ 10 M , leaving behind dim remnants of dwarf galaxies (Dekel & Silk 19861. Presumably 
the luminous dwarfs accrete gas at later epochs. Most gas is ejected by the first generations of SNe for 
systems with escape velocity ^ 50 km/s, leaving dim stellar remnants behind. 

In very low-mass halos gas cannot even fall in, because its specific entropy is too high ( |Rees| 



1986). This entropy barrier amounts to a temperature barrier since the gas density, which to first order 
is proportional to the total mass density, is the same in different halos at a given epoch. Only halos 
of mass £ 10 5 M trap baryons that are able to undergo early H 2 cooling and eventually form stars. 
Hydrodynamical simulations indicate that this lower limit is sharp ( |Gnedin 2000 ; Okam oto et al.||2008) . 
Reionization reinforces this limit by heating the intergalactic gas to high entropy, hence suppressing 
subsequent star formation (see Fig.[6]i. The abrupt increase of the sound speed to ~ 10 — 20 km/s at 
z ~ 10 means that dwarfs of halo mass ~ 10 6 — 10 7 M , which have not yet collapsed and fragmented 
into stars, will be disrupted. However massive dwarfs are unaffected, as are the high a peaks that develop 
into early collapsing, but rare, low mass dwarfs. 

2.5 Specific SFR 

Other serious, not unrelated, problems arise with low mass galaxies. In the hierarchical approach, these 
genetically form early. Theoretical models, both SAMs and hydrodynamical, appear to fail to account 
for the observed specific star formation rates (SFR per unit stellar mass or SSFR, Weinm ann et al.||2012| l, 
producing too little star formation at late times. Metallicity-dependent star formation alleviates the high 
redshift problem, reducing the stellar mass that is in place early and enhancing the SSFR as needed 
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Fig. 5 3D view of scaling relations of spheroidal systems from globular clusters (GC) to 
clusters of galaxies (CSph), via ultra-compact dwarfs (UCD), dwarf spheroidals (dSph), dwarf 
ellipticals (dE) and giant ellipticals (E), where the axes are half-luminosity, half-luminosity 



radius and total mass within half-luminosity radius ( |Tollerud et al. 201 1[ ). The red and blue 



planes respectively represent the Fundamental Plane and the "virial plane" of constant M/L. 



(Krumholz & Dekel) [2012 1. However, it leads to inconsistency at low redshift, because the change in 
metallicity and the gas fraction anti-correlate, hence leading to too little evolution in the SSFR. 

As shown in Fig.[7j the star formation time-scale (or 1/SSFR) goes from the MW value of ~ 10 Gyr 
at low redshift to ~ 0.5 Gyr at z ^ 2. This result suggests two distinct feedback-regulated modes of 
star formation: at low redshift via SNe and without AGN, and at high redshift with, most plausibly, 
quenching and possibly triggering by AGN playing a central role. One would expect a transition between 
these two modes as the AGN duty cycle becomes shorter beyond z ~ 1. 



A related triggering mechanism appeals to enhanced rate of merging at high z (Khochfar & Silk 



|201 1[ ). Alternatively, it has been argued that intensified halo cold gas accretion at early epochs may 
account for all but the most the extreme SFRs at high z, although this may require an implausibly high 
SFE HDekel et"aTl[2009l l. 



2.6 Spheroidal galaxies 

The baryon fraction is far from its primordial value in all systems other than massive galaxy clusters. 
SNe cannot eject significant amounts of gas from massive galaxies. Baryons continue to be accreted over 
a Hubble time and the stellar mass grows. One consequence is that massive galaxies are overproduced 
in the models, and that the massive galaxies are also too blue. 

Galaxies like the MW have peanut-shaped pseudobulges, in contrast with the classical bulges of 
more massive spirals. If formed by secular gas-rich disk instabilities, they should have an age distribution 
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Fig. 6 Evolution of circular velocities (at virial radius) for halos with efficient star formation 
(Ma mon et al.| |2012") . The smooth curves indicate the mean evolution of halos (with final 
masses log hM = 8 to 12 going upwards). The blue broken line is a model for the evolution 
of the minimum mass for galaxy formation (set by entropy feedback and related to the tem- 
perature of the IGM) and the cyan shaded region represents our ignorance of this parameter. 
The magenta curve is the maximum circular velocity for efficient gas infall. The grey shaded 
bands represent regions of thermal instability. The dashed curves shows 1 and 2 a fluctuations 
from the ACDM primordial density fluctuation spectrum. 



similar to that of the old disk. However the formation time would be at least ~ 1 Gyr. The elevated 
a/[Fe] ratio of our bulge favors a shorter formation time. This would be more consistent with an early 
disk instability phase reminiscent of that associated with clumpy gas-rich galaxies observed at z ~ 
2. Massive clump merging provides a possible solution for forming bulges at a relatively late epoch 
( |Ceverino et"aTl[20T0l ). 

However the time-scale (several Gyr) is too long to result in the enhanced a/[Fe] ratios character- 
istic of massive spheroids, or even the less extreme enhancement in the MW bulge. The shorter time- 
scales required arise in more plausible cosmological initial conditions that result in a redshift z > 2 for 



pseudobulge formation (Okamoto 2012 1. 



2.7 The role of AGN 



SNe have little impact on the formation of massive galaxies. Feedback from SN explosions fails to stop 



the streaming of cold flows towards the centre ( Powell et al. 201 1 1. The SN ejecta tends to be driven out 
with only modest interaction with, and entrainment of, cold infalling gas. A more coherent and effective 
interaction is provided by AGN feedback from supermassive black holes (SMBH). A clue towards a 



solution for these dilemmas comes from the accepted explanation of the Magorrian et al. relation, which 



relates SMBH mass to spheroid mass (Magorri an et al.| [1998) and velocity dispersion ( |Ferrarese &| 
Merritt||2DD0l see Fig. |g). 
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Fig. 7 Evolution of the specific SFR (SSFR) of galaxies of stellar mass 0.2 - 1 x 10 10 M 
(Weinmann et al. 201 1| top); SFR for galaxies from different samples (different color sym- 



bols), highlighting the "Main Sequence" and a population of starbursts, with SSFR mass- 
dependence inset ( Rodighiero etaL 2011 bottom). 



This requires collusion between black hole growth and the initial gas content of the galaxy when 
the old stellar spheroid formed. One conventionally appeals to outflows from the central black hole that 
deliver momentum to the protogalactic gas. When the black hole is sufficiently massive, the Eddington 
luminosity is high enough that residual gas is ejected. An estimate of the available momentum supply 
come from equating the Eddington momentum with self-gravity on circumgalactic gas shells, £Edd/c = 
AttGM/k — GMM gas /r 2 , where n us the opacity. Blowout occurs and star formation terminates when 
the SMBH-cr^ relation saturates. This occurs for Mbh cr^, close to the observed slope of > 5 
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Fig. 8 Black hole mass versus spheroid velocity dispersion (luminosity-weighted within one 
effective radius), from |McConnell et al.| | |20lT) 



(Graham et al. 201 and gives the correct normalization of the relation, at least in order of magnitude. 
This is the early feedback quasar mode. 

There is also a role for AGN feedback at late epochs, when the AGN radio mode drives jets and 
cocoons that heat halo gas, inhibit cooling, resolve the galaxy luminosity function bright end problem 
and account for the red colors of massive early-type galaxies. AGN feedback in the radio mode may 



also account for the suppression in numbers of intermediate mass and satellite galaxies (e.g., Cattaneo 



et al. 2009 and references therein). Feedback from AGN in the host galaxies also preheats the halo gas 
that otherwise would be captured by satellites. 



2.8 Galaxies downsize 



Our understanding of galaxy formation is driven by observations. Prior to 2000 or so, it was accepted 
that hierarchical galaxy formation predicted that small galaxies form prior to massive galaxies. The first 
indications that this was in error came from the recognition that more massive early-type galaxies have 
redder colors ( de Vaucouleurs 19611), higher metallicities fFaber 1 973 i and enhanced [a] / [Fe] metallic- 



ity ratios (Ziegler et al. 2005 1, indicative of an older stellar population with a shorter star formation time 
(see Fig.|9jl. This effect is called downsizing, as the most massive galaxies have their stellar populations 
in place early. In effect, we have a cosmic clock: incorporation into stars of debris from SNe II ( < 10 8 
yr) versus SNe I ( 10 9 yr) provides a means of dating the duration of star formation. This result was 
soon followed by infrared observations that showed that stellar mass assembly favored more massive 



systems at earlier epochs ( Gonzalez et al. 201 la I 
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Fig. 9 Metallicity ratio and age versus galaxy velocity dispersion (i.e. mass) and history of 
star formation (Thoma s et al.||20T0| l 



2.9 Morphological evolution 

We cannot do justice in this review to a largely phenomnological discussion of morphological evolu- 
tion of both disk and irregular galaxies. Here, observations are far ahead of theory. However, there are 
strong arguments to support a continuous sequence between dwarf spheroidal galaxies and SO galaxies 
(Kor mendy & Bender||2012] l. The transformation applies to the disk components and may involve ram 
pressure stripping of cold gas (Gunn & Gott 1972| l as well as galaxy harrassment (Mo ore et al.||1998| l. 



This sequence seems to acts in parallel to the pseudobulges or bulges of SO galaxies being generated via 
stripped/harassed or simply starved disk galaxies (Kormendy & Bender, 2 012| l. 



3 METHODS 

3.1 Observational surveys 

The fundamental driver of progress in astronomy is through observations. The advent of large galaxy 
surveys, either wide spectroscopic surveys probing the nearby Universe (e.g., SDSS) or narrower surveys 
using photometric redshifts and often in the infrared domain (e.g., with Spitzer and Herschel) to probe 
distant galaxies in the optical and near-infrared domains, has led to formidable progress in understanding 
galaxy formation. Nevertheless, it is difficult to link the galaxies we see at high redshift with the ones 
we see in local Universe, and one is prone to Malmquist bias, as well as aperture and other selection 
effects. 
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3.2 Semi- Analytical Models 

Several simulation techniques have been developed to be able to link galaxies from the past to the 
present, and to obtain a statistical view of the variety of the evolution histories of galaxies, in terms of 
star formation, stellar mass assembly and halo mass assembly. 




Fig. 10 Illustration of halo merger tree (Lacey & Cole 1993 i showing the progenitors of a 
halo selected at time to 



Given current computational constraints, it is impossible to achieve the sub-parsec or finer resolution 
needed to adequately model star formation and accretion onto black holes in a cosmological simulation. 
Theorists have invented a swindle, wherein the complex processes of star formation and accretion onto 
SMBH are hidden inside a black box called "sub-grid physics" that can be tagged onto a large-scale 
simulation. In SAMs, galaxies are "painted" on halos built from halo merger trees or detected in cos- 
mological dissipationless (dark matter only) simulations. The former (see Fig. 10 1 produce the mass 



assembly history (MAH) of halos with the condition that they end up in a halo of mass Mo at epoch zq 
(usually zq — 0). The branches are drawn from random samplings given the known conditional prob- 
abilities arising from extensions ( Bower] [l 99 l||Lacey & Colej |1993) > and modifications of the |Press"&| 
Schechter ( 1974\ formalism. Halo merger trees have the advantage of being rapid to compute, but lack 
positional information. Cosmological simulations, with up to 10 8 particles, are becoming increasingly 



common, but are intensive to process, in particular to detect halos (Knebe et al. 2011| and references 
therein) and subhalos ( |Onions et al.||2012| and references therein) and build halo merger trees. 

Once the halo MAH is identified, one follows the branches of the tree from past to present, to build 
galaxies. The galaxy formation recipe includes several ingredients: 

1. The gas cooling time must be short for the gas to dissipatively cool into a disk. In particular, gas 
cannot fall onto low-mass halos because of the cooling barrier and falls less efficiently onto high- 
mass halos because of a virial shock, whereas gas can infall along cold filaments on lower mass 
halos. 

2. Star formation occurs at a rate m = est m &as /t& yn , where t^ yn is a measure of the dynamical time 
of the galaxy. 

3. Feedback from SNe and from the relativistic jets arising from central SMBHs hiding as AGN heats 
up the surrounding gas. 
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4. While the gas settles into disks, major mergers of galaxies cause disks to transform into ellipticals, 
and after subsequent disk build-up, the merger remnant is identified to a bulge inside a spiral galaxy. 
The bulge can also be built-up by repeated minor mergers, as well as starbursts and secular evolution 
of the disk. 

5. The SAM keeps track of star formation times and predicts galaxy luminosities in different wave- 
bands using population synthesis codes. 

6. When a smaller halo enters a larger one, it becomes a subhalo, its galaxy becomes a satellite, and 
usually, the gas that attempts to fall onto the subhalo is now directed towards the central galaxy of 
the halo. 

7. Satellite galaxy trajectories are assumed to be those of the subhalos they belong to, and when they 
are no longer resolved in the cosmological simulation, or if one is only using a halo merger tree, 
the galaxies are merged with the central galaxy on a dynamical friction timescale (calibrated on 
simulations, e.g. |Jiang et al. 2008 ). 



An advantage of SAMs is that it is easy to gauge the importance of various physical processes by 
seeing how the outcome is changed when a process is turned off in the SAM. 

Present-day SAMs are increasingly complex, and SAMs can include up to 10 5 lines of code. The 
popularity of SAMs has increased with public-domain outputs (Bower et al. , 2006', Croton et al. 2006; 
|De Lucia & Blaizot| [20071 |Guo et al.||20"TT] ) and codes ( [Benson] |20i2] i~ 

3.3 Hydrodynamical simulations 

The weakness of SAMs is that much of the physics is controlled by hand (except for gravity, when 
the SAMs are directly applied to cosmological N-body simulations of the dark matter component). 
Hydrodynamical simulations provide the means to treat hydrodynamical processes in a much more self- 
consistent manner, and cosmological hydrodynamical simulations have been run for nearly 25 years 
(starting from |E ward 1988| l. 



These simulations come in two flavors: cell-based and particle-based. It was rapidly realized that 
cell-based methods could not resolve at the same time the very large cosmological scales and the 
small scales within galaxies, and the early progress in the field was driven by the Smooth Particle 
Hydrodynamics (SPH) method ( |Gingold & Monaghan[ [T9T7] |Monaghan| [19921 |Springel| |2010b| l, in 
which the diffuse gas is treated as a collection of particles, whereas the physical properties (tempera- 
ture, metal content, etc.) are smoothed over neighboring particles using a given SPH-smoothing ker- 
nel. Despite early successes (comparisons of different hydrodynamical codes by Frenk et al. 1999), 



SPH methods fail to resolve shocks as well as Rayleigh-Taylor and Kelvin-Helmholtz instabilities 



(Scannapieco et al. 2012[ ). This has brought renewed popularity to cell-based methods, with the ma 



jor improvement of resolution within the Adaptive Mesh Refinement (AMR) scheme (Kravtsov et al. 
|1997[ |Q'Shea et al.| |2005[ |Teyssier[ [2002 ), where cells can be refined into smaller cells following a 
condition on density or any other physical property. Moreover, schemes with deformable cells (that do 
not follow the Cartesian grid) were developed 20 years ago, and are now becoming more widely used 



(e.g., AREPO,|Springel 2010a). 



In these hydrodynamical codes, stars can be formed when the gas is sufficiently dense, with a 
convergent flow and a short cooling time. Current codes do not have sufficient mass resolution to resolve 
individual stars, so the star particles made from the gas are really collections of stars, with an initial mass 
function. One can therefore predict how many core-collapse SNe will explode after the star particle 
forms, and the very considerable SN energy is usually redeposited into the gas by adding velocity kicks 
to the neighboring gas particles and possibly also thermal energy. Similarly, AGN can be implemented, 
for example by forcing a Magorrian type of relation between the SMBH mass and the spheroidal mass 
of the galaxy, and feedback from AGN jets can be implemented in a similar fashion as is feedback 
from SNe. Hydrodynamical codes are therefore not fully self-consistent, as they include semi-analytical 
recipes for the formation of stars and the feedback from SNe and AGN. 
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3.4 New methods 

3.4.1 Analytical models 

The growing complexity of SAM codes (e.g., the publicly-available GALCTICUS code of Be nson] |2012| 
contains over 120000 lines of code and involves over 30 non-cosmological parameters with non-trivial 
values) has led some to seek simpler descriptions of galaxy formation. 

One level of simplicity is to parameterize the time derivatives of the different components of 



galaxies (stars, cold gas, hot gas, dark matter) as linear combinations of these parameters (Neistein 



& Weinmann 2010). But one can go to an even simpler level and characterize the fraction of gas that 
can cool (Blanchard et al.|[l992) or the mass in stars (Catt aneo et al.||2011| l as a function of halo mass 
and epoch. Although such approaches are much too simple to be able to capture the details of galaxy 
formation, they are sufficient to study simple questions. 

Assuming that all the gas in the range of temperatures between 10 4 K and the maximum where the 
cooling time is shorter than the age of the Universe (or the dynamical time) effectively cools, and that 
the gas is replenished one the timescale where halos grow, Blanchard et al7|( | 1992[ > showed that nearly all 
the baryons should have converted to stars by z = 0. Since this is not observed, this simple calculation 
shows that feedback mechanisms are required to prevent too high star formation. 




Fig. 11 Stellar versus host halo mass from analytical model by Cattaneo et al. (201 1 
dark matter simulation. The shaded regions are results obtained from the conditional luminos- 
ity function (blue, Yang et al. 2009) and abundance matching {gold, Guo et al. 2010). The 
large symbols denote galaxies who have acquired most of their stellar mass through mergers 
(rather than smooth gas accretion). The green curves show the galaxy formation model at 
z = and z = 3. 
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|Cattaneo et al.| ( |201 l| l apply their simple galaxy formation prescription onto the halos of a high- 
resolution cosmological N body simulation and reproduce the z = observed stellar mass function 
with only four parameters, despite the overly simplistic model. They find a fairly narrow stellar versus 
halo mass relation for the dominant ("central") galaxies in halos and a gap between their stellar masses 
and those of the satellites, in very good agreement with the relations obtained by |Yang et al.| (|2009 ) from 
the SDSS using conditional stellar mass functions (Fig. Hi. This gap is remarkable, as it is less built-in 



Cattaneo et al. s method than it is in SAMs: for example, halos with two-dominant galaxies (such as 
observed in the Coma cluster) are allowed. Similar "successful" analytical models have been proposed 



by Peng et al. (2010i and Bouche et al. (2010). 



3.4.2 Halo Occupation Distribution 

A simple way to statistically populate galaxies inside halos, called Halo Occupation Distribution (HOD) 
is to assume a functional form for some galaxy statistic in terms of the halo mass. The galaxy statistic 
can be the multiplicity function (for galaxies more massive or more luminous than some threshold, 

see upper left panel of Fig. 



Berlind & Weinberg 2002 



12i 



generically denoted CLF for conditional luminosity function (Yang et al 



the luminosity or stellar mass function, 
2003). Although these HOD 



methods have no underlying physics, they are a very useful tool to derive galaxy trends with halo mass, 
or, in other words to find the effects of the global environment on galaxies. 



3.4.3 Abundance Matching 

An offshoot of HOD models is to link the mean trend of some galaxy property in terms of the mass of 
its halo, using so-called Abundance Matching (AM). The idea is to solve N(> x) = N(> A/h a io)> i.e. 
matching cumulative distributions of the observed galaxy property, x, with the predicted one for halo 
masses, determined either from theory ( Press & Schechter |1974| Sheth et aL] 200 l|l o r from cosmo 



2006||Tinker et al.| |2008||Crocce et al.||2010||Courtin et al. 



logical N body simulations (Warren et a 
[20TTT ). 

Common uses of AM involve one-to-one correspondences between 1) stellar and halo mass for 
central galaxies in halos, 2) total stellar mass and halo mass in halos, and 3) stellar and subhalo mass in 
galaxies. Marinoni & Hudson (2002 1 performed the first such AM analysis to determine Al\ la \ /L versus 
L; they first had to determine the observed cosmic stellar mass function, not counting the galaxies 
within groups, but only the groups themselves (Marinoni et al. 2002). Guo et al. ( 2010| l used this 
third approach (called subhalo abundance matching or SHAM, and pioneered independently by | Vale &| 
Ostriker 2006 and Conroy et al. 2006 1, to determine the galaxy formation efficiency m stars /Afh a i as a 



function of Mh a i , by matching the observed stellar mass function with the subhalo mass function that 



they determined in the Millennium (Springel et al. 20051) and Millennium-II (|Boylan-Kolchin et al. 



2009) simulations. Although AM methods are based upon a fine relation between stellar and halo mass 
they can easily be adapted to finite dispersion in this relation (Behroozi et al. 2010 1. 

Not only is AM very useful to determine, without free parameters, the relation of stellar to halo 



mass (lower panel of Fig. 12 1, but it superbly predicts the galaxy correlation function of SDSS galaxies 
( |Conroy et al.| |2006| upper right panel of Fig. [T2| . The drawback of AM methods is that they do not 
clarify the underlying physics of galaxy formation. 

4 RESULTS FROM NUMERICAL SIMULATIONS 

4.1 General results from semi-analytical models of galaxy formation 

SAMs have been remarkably successful in constructing mock catalogs of galaxies at different epochs 
and are used in motivating and in interpreting the large surveys of galaxies. For example, they reproduce 
very well the z = stellar mass function and correlation function (see Fig. 13 ). 

However, attempts to solve the problems of high redshift galaxies have so far been woefully in- 
adequate. For example, they cannot reproduce the rapid decrease in the cosmic SFR since z = 1 (see 
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Fig. 12 Top left: Illustration of HOD models of multiplicity functions (per halo) obtained 
from abundance matching ( Conroy et al.| [2006) . Top right: Abundance matching prediction 
on the galaxy correlation function compared to SDSS observations (symbols), while the halo- 
halo correlation function is shown as dotted lines (Conroy et al. [2006). Bottom: Halo mass 
for given stellar mass obtained by abundance matching (AM), HOD, conditional luminosity 
function (CLF) and group catalogs (CL) at z = 0.1 ( |Behroozi et al.| [2012). The shaded region 
shows the AM analysis of |Behroozi et al.| ( |2010| l. 



Fig. 14 1. The early SAM feedback models used AGN quenching, and required excessive dust in early 



types in the nearby universe (Bower et al. 2006} . Refinements to high redshift attempted to account 
simultaneously for galaxy and AGN accounts, and only succeeded by requiring excessive amounts of 
dust in order to hide most of the AGNs seen in deep X-ray surveys (Fanidakis et al. 201 1 ). An early 
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Fig. 13 Illustrations of predictions of SAMs at z — 0. Upper left: Stellar mass functions (Guo 



et al.| |2011| l: symbols are from SDSS ( |Li & White] |2009| l, while curves are from a SAM 



run on both wide and low-resolution Millennium Simulation and on the higher but smaller 
MS-II simulation. Upper right: Galaxy correlation functions ( |Guo et aL]|2011| l. Bottom left: 
Evolution of the cosmic SFR (|Guo et al.| 201 l| l. Bottom right: Very low end of the galaxy 



luminosity function ( Koposov et al. 2009 1 



indication that SAMs were entering uncertain territory can be seen in the early predictions of the cosmic 
star formation history: as numerical resolution was increased, the predicted SFR increased without limit 
( |Springel & He rnquist 2003). This makes one begin to doubt the predictive power of SAMs. 

Clearly, baryon physics is far more complicated than assumed in the early SAMs of the 1990s. In 
fact, we still lack an adequate explanation for the evolution of the stellar mass function. Attempts to 
patch up the problem at low redshift, to avoid an excess of massive galaxies, exacerbate the inadequacy 
of the predicted numbers of massive galaxies at high redshift (Fontanot et al. 2009 1. One attempt to 




Fig. 14 Evolution of stellar mass functions predicted by Guo et al. (2011 



and red circles represent observations by Perez-Gonzalez et al. ( 2008 i and 



. Open triangles 



Marchesini et al. 



(2009), respectively. Black and green curves represent the predicted stellar mass functions of 
galaxies, respectively before and after convolving the stellar masses by 0.25 dex measurement 



errors. 



correct the problem at large redshift incorporates for the first time thermally pulsing AGB (or carbon) 



stars in the models, and the extra NIR luminosity reduces the inferred galaxy masses (Henriques et al. 



201 1 1. However the price is that the lower redshift galaxy count predictions no longer fit the models. 



4.2 Feedback and dwarfs 

Dwarf spheroidal galaxies are dark matter laboratories, dominated by dark matter. However the numbers 
defy interpretation. Feedback is readily adjusted to reduce the numbers of low mass dwarfs (Koposov 
et al.| |2009|l, b ut the most massive dwarfs predicted by ACDM simulations are not observed (Boylan- 
Kolchin et al.[ |2012| l. This may be a function of the neglect of baryons in the Aquarius simulations: 



inclusion of baryons reduces the central densities of massive dwarfs (Zolotov et al. 2012| l. Unorthodox 
feedback (AGN) may also be a solution (Boylan-Kolchin et al. 201 1 1. Moreover, most low-mass dwarfs 
have cores rather than the cusps predicted by CDM-only simulations. Baryonic feedback may reconcile 



data on dwarf core profiles with simulations that include star formation and gas cooling (Oh et al. 201 1 



Governato et al. 2012| l, who find that SN-driven outflows help flatten dark matter central density cusps 
As mentioned earlier, enhanced early star formation and SN production creates strong tensions with the 
need for strong late low mass galaxy evolution. SN feedback at later epochs may turn cusps into cores 



by sloshing of more recently accreted gas clouds (Mashchenko et al. 2006), more recently addressed 
i Pontzen & Governato| (2012 1, who consider bulk gas motions and require short intense bursts of star 



formation. There may be evidence for such phenomena in dwarf galaxies ( Weisz et al. 2012 1 

Multiphase simulations (Pow ell et al.||201 1) confirm the effectiveness of SN-driven winds, but find 
that they do not lead to baryon ejection. In a multi-phase medium with more realistic filamentary accre- 
tion, outflows are only typically < 10% of the gas accretion rate. It is not clear whether SN feedback 
may still provide enough momentum to yield an acceptable fit to the low mass end of the galaxy lu- 
minosity function for the classical dwarfs. Ram pressure stripping ( [Gunn & Gott| |1972| |Mayer et al.| 
2007[ ) remains an alternative or complimentary mechanism, and morphological transformation of disks 
into dwarf spheroidals may be accomplished by repeated rapid encounters, i.e. "harassment" ( |Moore 
et al.||1998] l or gravitationally-induced resonances ( jD'Onghia et al. 2009 ) . 

SN feedback enables present day disk galaxy properties to be reproduced, including the Tully- 
Fisher relation and sizes, except for massive disks. More energetic feedback, from an AGN phase, is 
envisaged as a possible solution (McCarth y et al.|[2012[ ). Many galaxies, including early types, have 
extended star formation histories. Minor mergers provide an adequate gas supply to account for these 
(Kaviraj et al. 2009) 1. However, hydrodynamical studies of the baryonic evolution and SFR in low mass 



18 



Joseph Silk 1 ' 2 ' 3 , Gary A. Mamon 1 



galaxies disagree about whether or not one can reproduce their observed properties, including dark 
matter cores and baryon fraction. Outflows may reproduce the observed cores ( |Governato et al. 201 2} 
if the SFE is high at early epochs, but such models fail to result in the strong evolution observed at low 
redshift (Weinmann et al.||2012| l. 

Tidal disruption also plays a role in disrupting satellites whose orbits intersect the disk or bulge. 
Dramatic discoveries due to deep imaging of nearby galaxies with very small, wide field of view, 



telescopes confirm the ubiquity of tidal tails that trace dwarf disruption in the remote past (Martinez 
Delgado et al. 2010). Simulations provide a convincing demonstration that we are seeing tidal disrup- 
tion in action ([Cooper et al. 2010). An independent confirmation of disruption in action comes from 
studies of the tidal tails around the outermost MW globular star clusters such as Pal 13. Gaps in the tails 
(Grillmair 2009) indicate the presence of dark satellites. Numerical simulations ( |Yoon et aL 2011[ ) find 
that high M/L satellites of mass ~ 10 7 M are required, again a prediction of the CDM model. 




rotation velocity (km/s) 



Fig. 15 Baryonic Tully-Fisher relation (Mamon & Silk, in prep.). Symbols are from HI mea- 
surements, where the velocity is the flat part of the rotation curve {green, from McGaugh 



|2012| ) or from line-widths (black, Gur ovich et al.||20T0} magenta, \Ha\l et aL]|201 1 ). The grey 
line is the naive ACDM (slope 3) prediction with no feedback, while the brown dashed line is 



the (slope 4) prediction from MOND. Note that the inclination of Ho II is uncertain (Gentile 

|eTaLl[20T2l ). 



At z — 0, it is possible that SN feedback at intermediate and low masses combines with entropy 
feedback from photoionization at low masses to conspire to give a linear baryonic Tully-Fisher relation 
(BTFR), as observed (see Fig. 15 I. This is an important issue as the normalization, slope and linearity 
of the BTFR have been used as evidence for Modified Newtonian Dynamics (MOND, Milgrom 1983 ) 
and against ACDM. Indeed, McGaugh (2011) has pointed out that the observations of baryonic mass 
(stars plus cold gas) as a function of the velocity of the flat part of the rotation curve is very well matched 
by the MOND prediction (with no free parameters). He argues that considerable fine-tuning is required 
to bring the naive ACDM (slope 3) prediction with no feedback to match the data. Our best-fit model 
matches the data equally well (with three free parameters), but the entropy feedback (photoionization) 
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implies that the relation should curve at low masses, except if one considers galaxies in which the bulk 
of the stars formed before the reionization epoch (z > 6). |Dutton] (|2012 ) also matched the BTFR data 
with a SAM. Moreover, two of the lowest mass galaxies i n the|M cGaugh sample have rotation velocities 
corrected for asymmetric drift (see Begum & Chengalur, 2004 ), after which the rotation curve of these 
galaxies is roughly linear with radius, and therefore depends on the last data point obtained with radio- 



observations, contradicting the fiat part of the rotation curve sought by McGaugh 




Fig. 16 Left: Galaxy luminosity function at z = 8 (Bradley et al.| 2012 1. Middle and right: 



evolution of the galaxy luminosity function ( |Bouwens et al. 2012[ l and of its faint-end slope 
( [Bradley et al7||2&12> . 



Intermediate-mass dwarfs are present at high redshift and have a steep luminosity function (Bradley 
|et al.|[20T2] see Fig. [To*]). They may contribute significantly to the reionization of the Universe. 

4.3 Gas accretion versus mergers 

Star formation seems to be too complex to be simply gravity-induced. Merging and AGN triggering are 
culprits for playing possible roles. What seems to be progressively clear is that there are two distinct 
modes of star formation. One mode occurs without any intervention from AGN and is characteristic of 
disk galaxies such as the MW, on a time-scale of order at least several galactic rotation times. Another 
mode is more intense, occurring on a relatively rapid time-scale, and involves the intervention of AGN, 
at least for quenching and possibly for enhancement or even triggering. 

The most important aspect of star formation is the role of the raw material, cold gas. There are two 
modes of gas accretion, which may be classified as cold flows/minor mergers and major mergers/cooling 
flows. The former provide supplies of cold gas along filaments, the latter a source of hot gas which may 
cool and feed star formation. 

The cold flows occur in filamentary streams that follow the cosmic web of large-scale structure (see 
Fig. [T7| i, and include minor mergers via the dwarf galaxies that similarly trace the web ( |Dekel et al.| 



2009 ). Theory suggests that, at low redshift, gas accretion by cold streams is important, and that the cold 
streams are invariably clumpy and essentially indistinguishable from minor mergers of gas-rich dwarfs. 
Major galaxy mergers account for the observed morphological distortions that are more common at high 
z, and generally lead to cloud agglomeration, angular momentum loss and cooling flows that feed star 
formation (Bou rnaud et al.|[2011b| ). 

Observationally, one finds that cold flows are rarely if at all observed. This is presumably because 



of the small covering factor of the filaments (Stewart et al. 2011 Faucher-Giguere & Keres 201 l| l. 



Indirect evidence in favor of cold accretion comes from studies of star formation in dwarfs. The best 



example may be the Carina dwarf where three distinct episodes of star formation are found (Tolstoy 



et al. 2009 [ l. However at high redshift, major mergers between galaxies are common. Indeed, Ultra 



Luminous Infrared Galaxie s (ULIRGs), whos e SFRs are huge, are invariably undergoing major, often 
multiple, gas-rich mergers (Borne et al. 2000 1 and dominate the cosmic SFR history at z £ 2, whereas 
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Fig. 17 Mass flux map of a M v = 10 12 M halo at z = 2.5 from a hydrodynamical simula- 
tion ( |Dekel et al.||2009) . The circle denotes the virial radius. 



normal star-forming galaxies predominate at low (z 2) redshift ( Le Borgne et al. 2009 1. This certainly 
favors the idea of massive spheroid formation by major mergers. 

U sing, their analy t ical m odel of galaxy formation on top of a high-resolution cosmological simu- 



lation, Cattaneo et al. (2011 1 show that only in massive galaxies (m star s > 10 ft M Q ) do galaxy 



mergers contribute to the bulk of the stellar mass growth (see also Guo & White 2008 who analyzed a 



simulation with 1 1 times worse mass resolution) and these mergers are mainly 'dry' (gas-poor). As one 
goes to lower stellar masses (down to their simulation's effective resolution limit of 10 10 ' 6 ft _1 M Q ) the 
role of mergers sharply diminishes, suggesting, by extrapolation, that mergers are, in general, unimpor- 
tant for the mass growth of both these intermediate-mass galaxies and low-mass galaxies, for which the 
bulk of the growth must be by gas accretion. Nevertheless, among those rare intermediate-mass galaxies 
built by mergers, the growth in mass is mostly in 'wet' (gas-rich) and minor mergers. In particular, the 
non-dominant cluster galaxies, known to be mostly dwarf ellipticals, are rarely built by mergers. 

The sudden dominance of major mergers at high galaxy masses is confirmed by trends with stellar 
mass of the colors, color gradients and elongations of SDSS galaxies (Tremonti et al. 2004| Bernardi 
|et al.||20f0l see also |van der Wei et aL]|2009[|Thomas et al.||2010|). At lowe r masses (and low redshift), 



minor mergers are required to account for sizes and masses (McLure et al. 2012 Lopez-Sanjuan et al. 
|20T2l >. 

Herschel observations of the Main Sequence of galaxy formation (SFR versus stellar mass) suggests 
that starbursts, commonly associated with major mergers, are displaced to higher mass and SFR, but 
only account for 10% of the SFR density at z = 2 (Rodighiero et al. 2011 1. However, this conclusion 
depends critically on the ~ 100 Myr timescale assumed for the starbursts. If the starbursts had shorter 
duration, say 20 Myr, given their effective observation time of ~ 1 Gyr, they would account for as much 
as 50% of the star formation at z ~ 2. It is difficult to gauge independent estimates of starburst age, but 
for example the UV continuum flattening observed at high z for luminous star-forming galaxies favors a 
younger starburst age (Gonzalez et al. 201 lt>| > as would possible SED corrections for nebular emission. 
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4.4 Initial stellar mass function 



The IMF of stars forming in galaxies is usually treated as universal in galaxy formation modelling. 
There has been a recent flurry of papers finding evidence for a systematic steepening, from the Chabrier 
(2003) to Salpeter ( 1955] ) IMFs, in massive early type galaxies. From a spectral absorption line analysis, 
a correlation of IMF steepening with enhanced velocity dispersion, [Mg/Fe] and sodium abundance is 



reported by (Conroy & van Dokkum 2012| l. A similar result is reported for stacked massive galaxy 
spectra (Ferreras et al. 20121. The modeling of the internal kinematics of early-type galaxies using in- 
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Fig. 18 Stellar mass-to-light ratio inferred from kinematical modeling (after subtracting off 
the contribution of the DM component), normalized to Salpeter ratio inferred from stellar 



populations versus stellar M/L from kinematical modeling, for six DM models (Cappellari 

|eTaTl|20T2l >. 



tegral field spectroscopy provides evidence for steeper IMFs (regardless of many plausible assumptions 
on the DM) in increasingly more massive galaxies (Capp eiTari et al.||2012| see Fig. [T8|). Lensing plus 



gas kinematics provides evidence for a Salpeter-like IMF in several massive ellipticals ( Dutton et al. 
2012b ]). Ther e may also be a correlation of a steeper IMF with the densest massive galaxies (jDutton 



et al.| |2012a| l. All of these studies report increasing M/L with increasing spheroid velocity dispersion 



and a/[Fe]. 

The possible degeneracy between IMF and DM fraction and shape is a concern because the DM 
profile steepens as a consequence of adiabatic contraction. While, Cappellari et al. (2012 1 tried a variety 
of DM models that do not significantly influence their result (since they only probed the region where 
dark matter accounts for at best 20% of the mass), only one study (jSonnenfeld et al. 2012|l so far has 
cleanly broken the degeneracy with the dark matter profile: By using a double Einstein ring, Sonnenfeld] 



et al. found a strong case for a Salpeter IMF. The adiabatic contraction of the DM is within the range 



found by |Gnedin et aL| (f2004 ). 

The implications of a steeper IMF in massive galaxies for galaxy formation models remain to be 
explored. The increased efficiency of star formation required at early epochs will certainly provide 
further tensions with the need to leave a substantial gas supply at late epochs for the observed late 
evolution observed for low mass galaxies, as discussed below. 
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4.5 Feedback and AGN 



Quenching of star formation has been largely motivated by the apparent success of SMBH feedback in 
reproducing the scaling and normalization of the black hole mass-spheroid velocity dispersion Mbh — 



a v ) relation, as first proposed by Silk & Rees ( 1998). SAMs indeed demonstrate that AGN feedback 
is able to quench star formation in massive elliptical galaxies ( |Croton et al.| |2006 Bower et al. 2006; 
Cattaneo et al.| [20061 ISomerville et al.| |2008|). One can reproduce the fairly sharp cut-off in the bright 



end of the galaxy luminosity function ( Bell et aL]|2003[[Panter et aT)|2007| ). These SAMs do not require 
"quasar mode" AGN feedback with Eddington luminosities. 





Fig. 19 Simulations (in 32 kpc box) of AGN feedback at 14 (left) and 22 Myr (right) after 



the onset of the jet, in edge-on (top) and face-on (bottom) views of log density (Gaibler et al. 

I20TT) 



High resolution hydrodynamical cosmological simulations indeed show that while cold streams 



initially feed the black hole, transferring angular momentum to produce central disks (Dubois et al. 



20 12a 1 that become gravitationally unstable and feed the compact bulge through migration of clumps 



(see Bournaud et al. 



(Dubois etal. 2012b) 



201 la 1, the cold flows are eventually interrupted by AGN -driven super-winds 



However the physics of driving SMBH outflows is still not well understood. One issue is that 
momentum-driven winds fail to account for the normalization of the AIbh — &v relation (Silk & Nusser, 
|2010[[Debuhr et ak 2012 ), with the shortfall being about a factor of 10. This momentum deficit can be 
supplied by radio jet-driven outflows ( Wagner & BicknellJ |201 1| ), which also account for the observed 



high velocities of entrained cold gas (Wagner et al. 2012 (. Alternative or complementary possibilities 



possibly more relevant to radio-quiet quasars, include positive feedback from outflow-triggered star for- 
mation ( Silk & Norman[ 2009[ |Gaibler et aL)|2011| see Fig. [19} and energy-driven outflows ( jFaucher- 



Giguere & Quataert 2012 1. Nearby AGN show dense molecular rings surrounding circumnuclear rings 
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of star formation ( |Sani et al7| |2012) >, reminiscent of the simulated triggering of star formation (Gaibler 
[eTaLl[20TT] i. 

SMBHs are generally found to correlate with bulges rather than with disks, pseudobulges or dark 
halos (Ho, 2007] |Kormendy et al.||20lTj|Kormendy & Bender) 201 1 1, although disk galaxies appear to 
follow a similar Mbh — &v relation, albeit with more scatter (Graham et al. 201 1 1. This would simplify 
formation mechanisms, suggesting that bulges and SMBH grow together, perhaps self-regulating each 
other. Massive black hole growth at early epochs seems to be (just) achievable by gas accretion. Large 



cosmological simulations (Di Matteo et al. 2012 see also Li et al. 2007, Sijacki et al. 2009;Khandai 



et al. 2012) 1 have shown that primordial massive BHs can grow by cold filamentary infall, and acquire 
masses of up to several billion solar masses by z = 6 in the most massive halos (M V - 1T ~ 10 12-13 M ). 

Insight into black hole growth is provided by looking for extreme deviations in the Mbh — &v 
relation. Massive black holes seem to be in place at high redshift before spheroids (Wang et al. 201 
This is also the case for a nearby starbust galaxy containing an AGN but without any matching spheroid 
or indeed massive stellar component ( |Reines & Deller||2012 l. On the other hand, SMGs seem to contain 
relatively low mass black holes for their stellar content (Alexander et al. 2008 ). 



5 FUTURE PROSPECTS IN OBSERVATIONS 



A clue as to the nature of a possible solution may come from the fact that quasars also reveal luminosity 
downsizing. This translates into downsizing of central SMBH mass. One might be able to connect the 
two phenomena if feedback from AGN was initially positive and also a strongly nonlinear function of 
SMBH mass. Predictions of positive feedback include circumnuclear rings on 10-100 pc scales in star- 
forming AGN. These should be resolvable with ALMA, via both molecular lines that probe pressurized 
molecular gas and FIR fine-structure lines that probe the interplay of intense FUV radiation fields with 
photodissociation regions (PDRs and XDRs). 

More conventionally, the evidence for AGN quenching of star formation seems strong. Superwinds 
driven by AGN are capable of depleting the reservoir of star-forming gas over relatively short time- 
scales. However, questions remain as to the relative roles of AGN winds, jet-driven bubbles, SNe, and 
radiation pressure especially from OB star clusters. No doubt, JWST, as well as 30+meter telescopes 
such as ELT will complement HST by producing spectacular IR images of star formation, AGN and 
outflows at work. Accretion of neutral gas will be studied at high sensitivity by SKA. Ultimately one 
needs a spectroscopic survey akin to SDSS at z = 1 — 2 and this will be provided by the Subaru Prime 
Focus Spectrograph with optical and NIR capability. The next decade should bring a vast increase in our 
phenomenological understanding of the basic processes at play in galaxy formation and evolution. 



6 FUTURE PROSPECTS IN ASTROPHYSICAL THEORY 

Theory lacks adequate resolution and physics. Of course these issues are intricately connected. One 
needs to tackle baryon physics and the associated possibilities for feedback. Today, state-of-the-art cos- 



mological simulations of the MW with gas and star formation, such as the ERIS simulation (Guedes 
|et al.|[20TT) , provide only w 100 pc resolution. Hence, in current simulations, the gas and star formation 
physics is included in an ad hoc way, because of the resolution limitation. For example, while stars are 
known to form in the dense cores — of density > 10 5 cm~ 3 — of Giant Molecular Clouds, the current 
hydrodynamical simulations adopt SF thresholds of typically 1 cm~ 3 and always <J 10 2 cm -3 . Sharp 



increases of the SF density threshold result in moving the SF regions outside of the nucleus ( Teyssier 



et al. 20101. However, in reality, it is the unresolved subgrid physics that determines the actual thresh- 
old, if one even exists. Mastery of the required subparsec-scale physics will take time, but there is no 
obvious reason why we cannot achieve this goal with orders of magnitude improvement in computing 
power. 

For the moment, phenomenology drives all modelling. This is true especially for local star forma- 
tion. A serious consequence is that physics honed on local star-forming regions, where one has high 
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resolution probes of star-forming clouds and of ongoing feedback, may not necessarily apply in the 
more extreme conditions of the early universe. 

One issue that arises frequently is whether the perceived challenges to ACDM justify a new theory 
of gravity. From MOND (Milgrom 1983] > onwards, there are any number of alternative theories that 
are designed to explain certain observations. However, none can explain the ensemble of observations 
any better than ACDM, nor do they rely on solid physical grounds. But to the extent that any unex- 
plained anomalies exist, these are invariably at no more than the 2 a level of significance. It seems that 
such "evidence" is not adequate motivation for abandoning Einstein-Newton gravity. Indeed, while it is 
overwhelmingly clear that there are many potential discrepancies with ACDM, we have certainly not 
developed the optimal ACDM theory of galaxy formation: the current models do not adequately include 
the baryons nor do we reliably understand star formation, let alone feedback. Other MOND-related is- 
sues are reviewed in Famaey & McGaugh (201 1), including challenges raised by the apparent emptiness 
of local voids and satellite phase space correlations. However, we regard these as more a matter of ab- 
sorbing the significance of ever deeper galaxy and 21 cm surveys, on the one hand (for example, deep 



blind HI surveys show that gas-rich galaxies are the least clustered of any galaxy population [Martin 



et al. 2012| l, and on the other hand, of questioning the details of hitherto inadequately modelled bary 
onic physics, as developed for example in |Zolotov et al.| ( [2012| l. Whether appeal to alternative gravity is 
justified by inadequate baryonic physics is a question of judgement at this point. Here is a summary of 
many of these failures: we cite some key reasons why ACDM does not yet provide a robust explanation 
of the observations: we list below several examples that represent challenges for theorists. 



1. Massive bulgeless galaxies with thin disks are reasonably common (Kormendy et al. 2010| l. 



Simulations invariably make thick disks and bulges. Indeed, the bulges are typically overly mas- 
sive relative to the disks for all galaxies other than SOs. Massive thin disks are especially hard to 
simulate unless very fine-tuned feedback is applied. A consensus is that the feedback prescriptions 



are far from unique (Scannapieco et al. 2012| l. One appealing solution involves SN feedback. This 
drives a galactic fountain that feeds the bulge. A wind is driven from the bulge where star formation 
is largely suppressed for sufficiently high feedback ( |Brook et al.||2012) . Another proposal includes 
radiation pressure from massive stars as well as SNe. The combined feedback helps expand the halo 
expansion, thereby limiting dynamical friction and bulge formation (Maccio et al. 2012| l. 

2. Dark matter cores are generally inferred in dwarf spheroidal galaxies, whereas ACDM theory pre- 
dicts a cusp, the NFW profile. Strong SN feedback can eject enough baryons from the innermost 
region to create a core (Governato et al. 2010| Pontzen & Governato 2012 1, but this requires high 
early SN feedback or a series of implausibly short bursts of star formation. 

3. The excessive predicted numbers of dwarf galaxies are one of the most cited problems with ACDM. 
The discrepancy amounts to two orders of magnitude. The issue of dwarf visibility is addressed 
by feedback that ejects most of the baryons and thereby renders the dwarfs invisible, at least in 
the optical bands. There are three commonly discussed mechanisms for dwarf galaxy feedback: 
reionization of the universe at early epochs, SNe, and (ram pressure and tidal) stripping. AGN- 
driven outflows via intermediate mass black holes provide another alternative to which relatively 



little attention has been paid ( Silk & Nusser 2010 1 



None of these have so far been demonstrated to provide definitive solutions. Reionization only 
works for the lowest mass dwarfs. The ultrafaint dwarfs in the MW may be fossils of these first 
galaxies (as checked by detailed models Koposo v et al.| |2009 ; Salvadori & Ferrari] |2009| |Bovill| 



& Ricotti 201 1). It i s argued that SN feedback solves the problem for the more massive dwarfs 



(Maccio et al. 2010| >. However, this conclusion is disputed by Boylan-Kolchin et al. (2011), who 



use the Aquarius simulations (Springel et al. , 2008 ) to predict more massive dwarfs in dark-matter- 
only simulations than are observed. These authors argue that the relatively massive dwarfs should 
form stars, and we see no counterparts of these systems, apart possibly from rare massive dwarfs 
such as the Magellanic Clouds. We have previously remarked that omission of baryonic physics 
biases the dark matter-only simulations to an overstatement of the problem by overpredicting dwarf 
central densities ( |Zolotov et al.||20T2) l. 
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6. 



The SFE in dwarfs is highly debated. Let us put aside the high SFE at early epochs that is required 
to obtain strong feedback in order to generate cores. For example, it is possible that intermediate 
mass black holes could be invoked to solve this problem and simultaneously generate the required 
low baryon fraction (Peirani et al. 2012 1. 

In order to obtain the required late epoch evolution (Weinmann et al. , 2012), one might appeal 
to a lower SFE in dwarfs, plausibly associated with low metallicities and hence low dust and H 2 
content. Models based on metallicity-regulated star formation can account for the numbers and 
radial distribution of the dwarfs by a decreasing SFE (Kravtsov 2010| >. This explanation is disputed 
by |Boylan-Kolchin et al.| ( |201 1| ), who infer a range in SFEs for the dwarfs of some two orders of 
magnitude. A similar result appeals to varying the halo mass threshold below which star formation 
must be suppressed to account for the dwarf luminosity function, whereas the stellar masses of many 
observed dwarfs violate this condition (Ferrero et al. 201 1 1. Finally, tidal stripping may provide a 
solution (Nick erson et al.[|201 1) , at least for the inner dwarfs. 

Another long-standing problem relates to downsizing. Massive galaxies are in place before lower 
mass galaxies as measured by stellar mass assembly, and their star formation time-scales and chem- 
ical evolution time-scales at their formation/assembly epoch are shorter. One popular explanation 
( Cattaneo et al.||2008[ ) is that galaxies cannot accrete/retain cold gas in massive halos, either be- 
cause of AGN feedback or because of virial shocks that prevent the gas supply of the disk in cold 
filaments ( |Birnboim & Dekel|[2003] l. 

It is possible to develop galaxy formation models with suitable degrees and modes of feedback 
that address many of these issues. However, a major difficulty confronted by all SAMs is that the 
evolution of the galaxy luminosity function contradicts the data, either at high or at low redshift. 
The SAMs that are normalized to low redshift and tuned to account for the properties of local 
galaxies fail at high redshift by generating too many red galaxies (Fontanot et al. 2009). Too few 
blue galaxies are predicted at z = 0.3. This problem has been addressed by including AGB stars 
in the stellar populations. This fix results in a more rapid reddening time-scale by speeding up the 
evolution of the rest-frame near-infrared galaxy luminosity function (Henriques et al. 2011| ). There 
is a price to be paid however: now there are excess numbers of blue galaxies predicted at z = 0.5. 
There is a well-known difficulty in matching both the galaxy luminosity function and Tully-Fisher 
scaling relation, even at z = 0. Reconciliation of the Tully-Fisher zero point with the galaxy lu- 
minosity function requires too high an efficiency of star formation (Guo et al. 2010). In fact, the 
problem is even worse: the models of massive spirals tuned to fit the Tully-Fisher relation are too 
concentrated ( |McCarthy et al. 2012| ). This is a reflection of the over- massive bulge problem in disk 
galaxies that simply refuses to go away ( |Navarro & Steinmetz 2000; Abadi et al. 2003 i. 
The luminosity function problem is most likely related to another unexplained property of high 
redshift galaxies. The SSFR evolution at high z is very different from that at low z. Essentially, 



10. 



11. 



it saturates. One finds an infrared Main Sequence of galactic SFRs: SFR versus M* (Elbaz et al. 



201 1 1. Neither the slope nor the scatter are adequately understood. Starburst galaxies lie above the 



Main Sequence, but the fraction of cosmic star formation in these systems depends on inadequately 
justified assumptions about starburst duration. For example, nebular emission and dust extinction 
affect infrerred ages, and one cannot easily understand the blue continuum slopes oberved at high 
redshift and lower UV luminosities (Bo uwens et al.||2~01 \\ . 

The observed rapid growth of early-type galaxy sizes since z = 2 for fixed stellar mass cannot be 
reproduced in SAMs or analytical models (Cimatti et al. 2012| ): at z — 2 galaxies are too compact. 
Much has been made of nearby rotation curve wiggles that trace similar dips in the stellar surface 
density that seemingly reduce the significance of any dark matter contribution. Maximum disks 
optimize the contribution of stars to the rotation curve, and these wiggles are most likely associated 
with spiral density waves. A similar result may be true for low surface brightness gas-rich dwarf 
galaxies ( Swaters et al. 201 1] >. 



High mass-to-light ratios are sometimes required for maximum disk models of spiral galaxy rotation 
curves, but these are easily accommodated if the IMF is somewhat bottom-heavy. The case for 
IMF variations has been made for several data sets, primarily for early-type galaxies (e.g., see 
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12. 



13. 



14. 



|Dokku m & Co nroy||2011| l. The LSB dwarfs are plausible relics of the building blocks expected in 
hierarchical formation theories. 

Spiral arms are seen in the HI distribution in the outer regions of some disks. This tells us that 
significant angular momentum transfer is helping feed the optical inner disk. The baryon self-gravity 
is large enough that one does not for example need to appeal to a flattened halo, which might 
otherwise be problematic for the DM model ( jBertin & Amorisco 2010| l. 

The slope and normalization of the baryon Tully-Fisher relation do not agree with the simplest 
ACDM prediction. The observed slope is approximately 4, similar to what is found for MOND 



(Mil grom| |1983[), whereas ACDM (without feedback) gives a slope of 3 ([McGaughj |20 1 1 1 12012|>, 



but fails to account for the observed dispersion and curvature. 

The baryon fraction in galaxies is some 50% of the primordial value predicted by light ele- 
ment nucleosynthesis. These baryons are not in hot gaseous halos (Anderson & Bregman 2010| l. 



Convergence to the universal value on cluster scales is controversial: convergence to the WMAP 



value is seen for X-ray clusters above a temperature of 5 keV ( Dai et al. 



20101, but could be as 



large as 30% even for massive clusters ( Andreon 2010 Scannapieco et al. 



20121. If the latter dis- 



crepancy were to be confirmed, one would need significant bias of baryons relative to dark matter, 
presumably due to feedback, on unprecedentedly large scales. 

15. The distribution of the MW satellite galaxies in a great circle (Lynden-Bell 1982) > is unexpected in 
the ACDM co ntext (|Kroupa et al. 2005) >. However, infall onto halos is not spherically symmetric 
(Aube rFet al.| |2004| l, and subhalos tend to lie in a plane (Libeskin d et al.| |2005|l. The details of 
the thickness of this plane remained to be settled (e.g., |Kroupa et al.||2010| versus|Libeskind et al. 
|20TTT >. 

16. There is a significant lack of galaxies in comparison with standard expectations in the Local Void 
close to the Local Group (Peebles 2007; Tikhonov & Klypin 2009). But it is not yet clear whether 
this region fairly low galactic latitude region has been surveyed as closely as other regions. 

17. Bulk flows are found over 100 Mpc scales that are about two standard deviations larger than 

201Q) >. The technique primarily uses Tully-Fisher and 



expected in ACDM (Feldman et al. 



Fundamental Plane galaxy calibrators of the distance scale. An X-ray approach, calibrating via the 



kinetic Sunyaev & Zeldovich 



( 1972] > effect (kSZE), claims the existence of a bulk flow out to 800 



Mpc (Kashlinsky et al. 2010 1. However the discrepancies with ACDM are controversial because 
of possible systematics. A recent detection of kSZE confirms pairwise bulk flows of clusters at 4 a 
and is consistent with ACDM pfand et al7]|2012| >. 



Several of these issues may be linked. For example, the analysis of Cappellari et al. ( 2012| > that the 
IMF is non-universal, with shallower (top-heavy) IMFs for galaxies of lower velocity dispersion, can be 
linked with the known relations between velocity dispersion and metallicity (e.g., |Allanson et aT 2009) 



to produce a relation between IMF and metallicity, which goes in the right direction: low-metallicity 
systems have top-heavy IMFs. Until now, observers assumed a universal IMF when deriving stellar 
masses. They have therefore overestimated the stellar masses of low-metallicity systems. We would like 
to think that this overestimation of A/* might explain at the same time the evolution of the cosmic SSFR 
and that of galaxy sizes. Indeed, at high redshift, galaxies are expected to be more metal-poor, and the 
overestimate of their typical stellar masses will lead to an underestimate of their SSFRs, relative to 
those of lower-redshift galaxies. Therefore, the cosmic SSFR may not saturate at high redshift, which 
will make it easier to fit to models. At the same time, if high redshift galaxies have lower stellar masses 
than inferred from a universal IMF, then for a given stellar mass, they have larger sizes than inferred, and 
the too rapid evolution of galaxy sizes (relative to models) might disappear. We propose that observers 
replace stellar mass by if-band rest-frame luminosity, which, if properly measured, can serve as a useful 
proxy for stellar mass, independently of any assumed IMF. 

In summary, it is clear that many problems await refinements in theoretical understanding. No doubt, 
these will come about eventually as numerical simulations of galaxy formation are refined to tackle 
subparsec scales. 
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